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Abstract 

We consider a last-passage directed percolation model in Z+, with i.i.d. weights whose 
common distribution has a finite (2 + p)th moment. We study the fluctuations of the 
passage time from the origin to the point (n, n '-"-'). We show that, for suitable a (de- 
pending on p), this quantity, appropriately scaled, converges in distribution as n — » oo to 
the Tracy- Widom distribution, irrespective of the underlying weight distribution. The 
argument uses a coupling to a Brownian directed percolation problem and the strong 
approximation of Komlos, Major and Tusnady. 

1 Introduction 

The concept of universality class plays a key role in statistical mechanics, making it possible 
to classify a huge variety of models and phenomena by means of well chosen scaling exponents. 
For instance, many growth models are expected to share similar properties which fall into the 
framework of the KPZ universality class - see for example the survey by Krug and Spohn 

EH 

In this note, we focus on the particular example of directed last-passage percolation. Let 
ojf\ i > 0, r > 1 be i.i.d. random variables. We consider directed paths in the lattice Zl, 
each step of which increases one of the coordinates by 1. For n > 0, k > 1, the (last-) pas sage 
time to the point (n, k) is defined by 

T(n,k)= max { V J r) \ , (1.1) 

where IT(n, k) is the set of directed paths from (0, 1) to (n, k). More precisely, 

II(n, k) = | {z\,Zz, ■ ■ .,z n+k ) e {l? + ) n+k : z\ = (0, 1), z n+k = (n,k), 

z J+1 - zj G {(0, 1), (1,0)} for 1 < j < n + k - 1 I. 
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When the underlying weight distribution is exponential or geometric, scaling exponents 
for this model are rigorously known. The deviation from a straight line of the optimal path 
to the point (n, n) is of the order n 2 / 3 (corresponding to the exponent £ = 2/3) j3J 
while the fluctuations of the passage time T(n, n) are of the order n 1 / 3 (corresponding to the 
exponent \ = 1/3)- In fact, one can give much more precise information: for the exponential 
distribution with mean 1, say, it is shown in [H] that 

n~ 1/3 [T(n,n) -An] -> F TW , (1.2) 

where -Ftw is the "Tracy- Widom" distribution, which also appears as the asymptotic distri- 
bution of the largest eigenvalue of a GUE random matrix. It is expected (but not yet proved) 
that the same scaling exponents, and indeed the same asymptotic distribution in (|1.2|) . should 
hold for a wide class of underlying weight distributions. 

In this note we give a universality result for the quantities T (n, \_n a \) f° r a < 1. Thus, we 
are concerned with passage times to points which are asymptotically rather close to the hor- 
izontal axis, for a general class of underlying weight distribution. Our result is the following: 



Theorem 1 Suppose that E |w} r ''| p for some p > 2, with li = Ewf, cr 2 = Var(w i (r) ). Then 
for alia < I 

T (n, [n a \ )-nfi- 2a rT^r ^ ^ 
an 2 ~ e 

in distribution. In particular, if the weight distribution has finite moments of all orders, then 
fO|) holds for all a < 3/7. 

Heuristically the theorem can be understood as follows. As the optimal path goes from 
the origin to (n, nL a J), one can imagine that between each step upwards, the path typically 
takes on the order of n 1 ~ a steps to the right. Thus it should behave as the optimal path from 
the origin to (n a , n a ) in a last percolation model with Gaussian weights of variance n 1_a . On 
the renormalized scale the expected fluctuations are of order (n ) 1 / 3 . In this way, we recover 
the fluctuation exponent 

I — a a 1 a 
X ~ 2 + 3 ~~ 2 ~ 6 ' 

This heuristic is made precise by coupling the discrete model with a Brownian directed 
percolation model for which the fluctuations have been explicitly computed. This is done 
using the strong approximation of a random walk by a Brownian motion due to Komlos, 
Major and Tusnady. (We note that this strong approximation has already been applied to 
similar last-passage percolation models by Glynn and Whitt 6 ). 

A different sort of universality result for paths near the axis in directed percolation models 
is given in |13) . By subadditivity, one has the convergence n~ 1 T(n, [a^J) — * j(x) a.s. and 
in C%, for some function 7. Under the hypothesis of Theorem^ it's shown that 7(3;) = 
u + 2a^fx + o(y/x) as x I 0. 

The proof of Theorem^is given in the next section. In Section|31 we make some comments 
on related models and possible extensions. 



2 



2 Fluctuations of the passage-time 

(r) 

We first introduce the Brownian directed percolation model. Let B\ , t > 0, r > 1 be a 
sequence of independent standard Brownian motions. For t > 0, fc > 1, define 

W(t,fc) = {(u ,wi,...,u fe ) € K fe+1 : = k <«i < ••• <u fc = t}, 

and then let 

k 

L( t) k)= sup ^ [S« -B^J. (2.1) 



ueu(t,k) r=1 

One can rewrite the definition of T(n, fc) at in an analogous way: 



T(n,k)= sup EteVi-^Lj]' ( 2 ' 2 ) 



u&J(n,k) r _j 

where = £™ q 1 w i • 

The random variable L(l, fc) has the same distribution as the largest eigenvalue of a fc x fc 
GUE random matrix @], 0, Hence in particular (e.g. |16|) 

fc 1/6 [L(l,fc)-2Vfc] ->F TW 

in distribution, where -Ftw is the Tracy- Widom distribution. 

By Brownian scaling, L(t,k) has the same distribution as i/tL(l,k). Using this we get, 
for any < a < 1, 

L ^J) ; 2^ fTw (2 . 3) 

712 6 

in distribution. 

Theorem^says that the same distributional limit as in 1)2. 3JI (in particular, with the same 
order of fluctuations) occurs for the law of T(n, [n a \), f° r & general underlying distribution 

(r) 

of the weights uj\ , if a is sufficiently small. We will use the following strong approxima- 
tion result, which combines Theorem 2 of Major an d Theorem 4 of Komlos, Major and 
Tusnady US]: 

Proposition 2 Suppose U3%, i = 1,2,... are i.i.d. with E \oJi\ p < oo for some p > 2, and with 
Kuii = 0, Var(aj.;) = 1. Let S m = J^Lq Wj, m > 1. 

Then there is a constant C such that for all n > 0, there is a coupling of the distribution 
of (ll>i, . . . ,0J n ) and a standard Brownian motion B t , < t < n + 1 such that, for all x € 
[nVp,n 1/2 ], 

P( max \B m - S m \ > x] < Cnx~ p . (2.4) 
ym=i,2,...,n+l / 

Proof of Theorem^} 

We may assume that /i = and cr 2 = 1, so that we need to prove that 

TMn-J)-^^ (25) 

n2 e 

in distribution (for general ^ and cr 2 , one can obtain l|1.3(l from 1)2. 5|l after replacing uj by 
(w-//)/cr). 
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If K (r) ) lir and (B[ r) ) tr are all defined on the same probability space, then from l|2.IJ) and 
we get 



\T(n, [n a \)-L(n, [n a ])\ 



sup 



Ln"J 

sup 

u'eW(n,|n a J) ; =1 



/ i \ u' r u' 



L» a J 



< sup 

ueW(n,L« a J) 



EW2 



(r) 



(r) _ B (r) 



5' 



(r) 
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(r) 



LB 



L«rJ+l 



B 



W _ o(r) 



(2.6) 



L™ a J 

r— 1 

L™ a J 



max 

,2,...,n+l 



S^-B^ 


+ sup 


B (r) _ B (r) 


> 




0<s,t<n+l 








|«— *|<2 




> 



r=l 

where we have defined 



r) 



max 

i=l,2,...,n+l 



and W^p = sup 

0<s,t<n+l 
Is— 1\<1 



For each r = 1, . . . , |/i a J we will couple ^Wq'"-' , w{ r \ . . . , Un'j and £j r \ 0<i<n + lasin 

Proposition |3 maintaining the independence for different r, so that the V« , 1 < r < \n a \ 
are i.i.d. with 

P (V„ (r) > a?) < Cnx-v (2.7) 

for all x S [n^n 1 / 2 ]. 

Let Ai be the event |maxi< r < [ n aj Vn^ > n 1 ^ 2 ^. Then from (|2.7|) . 



J (Ai) < rx°Cn(n 1 ' /2 )- p = C*n a+1 - p / 2 -> as n -> oo, 



(2.8) 



since by assumption a < | f | — |j < -P (l — p) = 1 — 1- 

Also let Ai be the event |max 1 < I ,<^ l aj Wn^ > Using the reflection principle and 

standard estimates on the normal distribution, 



\A 2 ) < n a t | sup 

0<s,t<n+l 
\s-t\<2 



> n 



l/p 



<«°VP sup B t - inf B t > n 1/p 

~^ \i<t<i+3 i<t<i+3 



< n a+1 P ( sup | B t | > n 1/p /2 

\0<t<3 
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= 4n a+1 p(B 3 > n 1/p /2\ 
< Cl n a+1 exp (-c 2 n 2/p ) 
— * as n — > oo. 

From H2.6fl . (|2.7[) and the definitions of the events A\ and A2, we have 



E 



|T(n, [n a \)-L(n, [n a \)\ ; Af U A c 2 

< 2n a E (V„ (1) + W { n 1] ; A% U A% 

< 2n a 



< 2n a Un 1 ^ + J P (V n W > a?) dx) 



< 2n a n 1/p + / Cnx- p dx 



= 2n a (n^ + Cznl-x-P+Xvl) 
< C 3 n a n 1/p , 

where C2 and C3 are constants independent of n. 

Together with 1)2. 8j> and 1)2. this gives that, for any e > 0, 



(2.9) 



(2.10) 



\T(n, \ n a \ ) — L(n, [n a \ ) | > n a+1/p+e j ^Oasn^oo. (2.11) 

The assumption a < f ^5 ~ jjjj implies that, for e sufficiently small, a+i + e < 5 — f • Thus 
|T(n, Ln°J)-L(n, L"- a J) 



in distribution, as 



n 2 b 

Using (|2.3[1 we obtain (|2.5|l as desired. 



(2.12) 
□ 



3 Further remarks 
3.1 Larger values of a 

It seems unlikely that the value 3/7 in Theorem Q] represents a real threshold. For a > 3/7, 
consider the typical difference between the weight of the maximal Brownian path and the 
weight of the discrete approximation; this will be large compared to the order of fluctuations 
of the maximum. However, the standard deviation of this difference may be smaller; one 
might expect it to be of order ?W 2 rather than order n a , since it is composed of n a terms of 
constant order which one expects to become independent as n becomes large. An argument 
along these lines would effectively allow us to replace n a by n a / 2 in (|2.1(J|) leading to a bound 
a < 3/4 rather than a < 3/7. However, above a = 3/4 it seems that the behaviour is genuinely 
different and more sophisticated arguments would be required: the fluctuations of the error in 
the discrete approximation to the Brownian path are likely to be larger than the fluctuations 
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of the maximal weight itself, and so one might no longer expect the maximal discrete path 
to follow closely the maximal Brownian path (even when the weights are "strongly coupled" 
to the Brownian motions as above). 



3.2 Transverse fluctuations 

The exponent x = \ — % which we obtain should be related to the transversal fluctuations of 
the optimal path away from the straight line {y = rt a_1 a;}. Let us introduce the exponent 



£ = lim 

n— >oc 2 logn 




where, say, Vi is the smallest value r such that the point (i,r) is contained in the optimal 
path. 

Corresponding to the universal value of \ obtained, one would expect that £ should be 
also a universal exponent equal to 2a/3. To see this, we follow the heuristics explained in 
the introduction. On a renormalized level, the optimal path should behave as the optimal 
path from the origin to (n a ,n a ) in a last percolation model with Gaussian weights. This 
would imply that the transverse fluctuations should scale like (n a ) 2 / 3 , where £ = 2/3 is the 
(predicted) standard fluctuation exponent for directed last-passage percolation. 

The strategy used in |H] for the derivation of the transverse fluctuations requires not only 
the knowledge of the last passage time fluctuations, but also a precise control of the moderate 
deviations. Our approach does not allow us to derive such sharp estimates (in particular 
we are missing some uniformity with respect to the direction of the path). For this reason, 
we do not yet have a proof of the universality of the transversal fluctuation exponent in our 
framework. 



3.3 Related models 

The last-passage percolation processes have a natural interpretation in terms of systems of 
queues in tandem (see for example |H E] )• Considering paths near the axis corresponds to 
considering regimes of very high or very low load in the queueing systems. In the case of 
exponential weight distribution, these queueing systems correspond closely to totally asym- 
metric exclusion processes or totally asymmetric zero-range processes. There are also close 
links with systems of non-colliding particles. See for example |14| for a survey. 

Of course, there are also strong connections between these models and random matrix 
theory. We mention one particular direction related to the topic of this paper. Let A n ^ be 
an n x k random matrix with i.i.d. entries, and let Y n ,u — A n _k{A n _k)* ■ In the special case of 
the Laguerre ensemble, where the common distribution of the entries is complex Gaussian, one 
has an explicit correspondence between the largest eigenvalue of Y and the passage time to 
(n, k) in a directed percolation model with exponential weights (see for example Proposition 

1.4 of [Hj, and [HJ and Section 6.1 of [2] for extensions). For a general distribution, there may 
not exist an explicit mapping between the matrix model and the directed percolation model, 
but we believe that a similar averaging mechanism to that observed in our context will also 



G 



play a role in the random matrix setting. Thus on the basis of the analysis of the last passage 
time fluctuations for paths close to the axis, we conjecture that when k and n tend to infinity 
with suitable rates, the fluctuations of the largest eigenvalue of Y n ,k should depend only on 
the mean and the variance of the coefficients of A n ^. 
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